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Abstract 

An analysis is presented that illustrates how the James Webb Space Telescope (JWST) fine-phasing 
process can be carried out using the Near-Infrared Spectrograph (NIRSpec) data collected at the science 
focal plane. The analysis considers a multi-plane diffraction model which properly accounts for the 
microshutter diffractive element placed at the first relay position of the spectrograph. Wavefront sensing 
results are presented based on data collected from the NASA Goddard Microshutter Testbed. 
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1. Introduction 

The JWST 1 (James Webb Space Telescope) is one of NASA’s great observatories of the 
Origins program ' and is scheduled to launch at the beginning of the next decade. The 
JWST will incorporate a multi-object spectrograph, or NIRSpec, 3 the Near Infrared 
Spectrograph, as one of the primary science instruments on-board the spacecraft. As 
such, the NIRSpec instrument will utilize a portion of the science focal plane, with /. in 
the range: 0.6 to 5 /x. to carry out investigations on galaxy formation, clustering, 
molecular structure, star formation, active galactic nuclei, and young stellar clusters. 

The JWST will be constructed using advanced optical technology and extremely light- 
weight mirrors. 4 The primary mir ror is a 6.5 meter diameter 18 segment hexagonal array. 
The figure quality should be diffraction limited at 2 /x. operate at 30-60 K. and have an 
areal density of less than 15 kg/m". As a result, the JWST commissioning phase and 
periodic maintenance/optical correction will utilize image-based wavefront sensing and 
control technology to align the mirror segments and minimize figure error as well as 
position the secondary mirror of the 3-mirror anastigmat design. The WFS method 
specified for JWST is image-based in the sense that point source images (or other known 
object) are collected to recover optical phase information. The primary camera for this 
function is NIRCam.' An important risk mitigation strategy for the JWST design is 
based on engineering redundancy for observatory functions. Therefore, an obvious 
redundancy to consider is the WFS process itself which logically implies the use of 
multiple cameras for WFS. In addition, it has been noted that imaging performance over 
the JWST FOV can be better compensated by using simultaneous WFS results from 
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multiple cameras using a Kalman filter approach. 6 Degradation in optical quality due to 
WFS and control at a single field point has been considered in. 7 

To address the concern for sufficient WFS redundancy as well as WFS and control 
induced anisoplanatism - one possibility is to augment NIRCam WFS data using focal 
plane images collected by NIRSpec (while operating in its camera imaging mode). 
However, NIRSpec WFS is complicated by the fact that the spectrograph design employs 
a microshutter 8 array at the first relay position of the spectrograph. As a result, diversity 
defocus data collected at the NIRSpec focal plane will be imprinted with a rectangular 
diffraction pattern induced by the microshutter array light shield. The subject of this 
paper is a discussion of WFS results obtained by processing these diversity defocus 
images using an image-based WFS algorithm. Modeling comparisons are also made 
utilizing multi-plane diffraction techniques which motivates the use of a simple overlay 
of the MS light shield upon the defocused intensity data. In Section 2, the microshutter 
arrays and GSFC Microshutter Testbed are discussed. Section 3 describes a multi-plane 
forward diffraction model that is used to calculate images for comparison to the MTB 
data. Image-based WFS is discussed in Section 4 the WFS results are presented. 
Conclusions are discussed in Section 5. 

2. Microshutter Arrays and Testbed Optical Design 

The microshutter array 8 is a field programmable micro-electromechanical device that 
enables the NIRSpec instrument to select hundreds of different objects in a single field of 
view. Microshutters are an alternative to the micro-mirror 9 systems which cannot meet 
JWST requirements for image quality, spectral coverage or contrast. Figure 1 shows a 
scanning microscope image of a portion of the microshutter array under development at 
the NASA Goddard Space Flight Center. 



Figure 1 . Scanning microscope image of a portion of the microshutters array. 

A microshutter array consists of about 1000 x 500 microshutters, 100 pi center-to-center 
and about 90 pi wide. NIRSpec will employ a 100 pi x 200 pi (spectral x spatial direction) 
array. All shutters are slightly magnetized and are opened by scanning a permanent 
magnet over the array. Selected apertures are held open electrostatically by applying a 
voltage difference between the shutter and an electrode on the wall. Following the single 
passing of the permanent magnet, the resilience of the hinges flips the remaining shutters 
closed. Additional reffences and information on the microshutter arrays can be found at 
http://bram.gsfc.nasa.gov/ms_webpage/ms_main.html. 



The MTB is a resource located at the GSFC and was designed to investigate and test the 
electro-mechanical operation of the shutters as well as electro-optical characteristics such 
as transmission and light leakage around closed shutters at infrared wavelengths. Early 
modeling efforts' 0 have focused on the electro-optical properties of the arrays. 


The MTB has the following layout as labeled in Figure 2: a tungsten point-source is 
filtered at 0.633 /x (and 2 /x) and then collimated by an off-axis parabola. The beam is 
then stopped by an iris which controls the plate scale at the microshutter array. This is 
variable to F/15 (the fastest). Lens assembly one is an IR/visible triplet which focuses the 
beam to the microshutter. CaF2 windows are used as windows to a dewar which houses 
the microshutter array. The beam is then re-collimated to a 2 nd stop and then focused 
using a 2 nd off-axis parabola. The testbed setup includes both an Apogee CCD with 13 /x 
pixels as and an Indigo Merlin 1.5 - 5.0 ft waveband InSb infrared Camera. The InSb 
array size is 320 x 256 with a 30 g pixel size. 
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Figure 2. MTB Optical Design. 

For comparison with the MTB layout. Figure 3 show's a MOS design in paraxial form. 
The point source at the MTB is analogous the input beam of the MOS design shown to 
the far left of Figure 3. 
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Figure 3. MOS in Paraxial Form. 

The off-axis parabola of the MTB collimates to pupil space and then lens assembly 1 
focuses the beam to the microshutter (labeled as (1) in Figure 3). The beam is re- 
collimated using lens assembly 1 (labeled as (2) in Figure 3) and then the final off-axis 
parabola of the MTB serves as the camera lens labeled as (4). As a result, these is direct 



correspondence between the MTB and a MOS design and therefore we expect some 
correspondence in the WFS data. For the WFS data, the IRIS-1 was at about 0.6” to 
make the beam approximately F/48 at the microshutter. This is then relayed to F/51 at 
the detector to give a slightly oversampled PSF by a factor of 1.25. 

3. Multi-Plane Diffraction Model and MTB Data 

Focus-diverse WFS requires intensity data collected at a defocused image plane. For 
NIRSpec, there are two options for introducing diversity defocus. The advantages and 
disadvantages to each are summarized in Table 1. 


Table 1. Options for implementing diversity defocus 


Diversity 

Defocus 

Secondary Mirror Motion 

NIRSpec Cam 

Advantage 

(1) Common path diversity 

(2) Large defocus possible 

(1) Secondary Mirror 
Stationary 

Disadvantage 

(1) Sec Mirror Motion; 

(2) Diffractive shutter or 
slit 

(1) too small de focus ? 

(2) camera movement 

(3) non-common path 
diversity 


The first option involves translating the secondary mirror along its optical axis. This has 
an important advantage in that all instruments of the ISIM science module will “see” 
diversity defocus and therefore WFS data can be collected in multiple cameras 
simultaneously. As a result, defocus consistency is easier to achieve since all cameras 
will receive common path defocus diversity. In addition, larger defocus values may be 
realized quite easily using secondary mirror motion. A disadvantage though is that 
defocus is introduced ahead of the MSA which has the consequence of introducing a 
rectangular grid structure onto the defocused image. However, as we show in this 
investigation, this artifact does not seriously degrade WFS performance. Another 
disadvantage is associated with the risk involved in moving such a major optical element 
of the OTE. Should the secondary mirror become stuck or frozen in its defocused 
position, the observatory will likely be rendered in-operable. 

The second option involves translating the NIRSpec imaging camera directly along the 
chief-ray optical axis. But with a limited range of motion it may not be possible to obtain 
a sufficient amount of de focus to perform accurate WFS. Diversity selection for image- 
based WFS has been previously analyzed based on the spatial frequency content of the 
input aberration. 11 Additional analysis on accuracy vs defocus has also been presented. 1 ' 
In addition, the problem of non-common path defocus diversity is introduced which 
would tend to aggravate the WFS&C induced anisoplanatism - which is one of the very 
reasons for using NIRSpec for WFS to begin with. However, translating the NIRCam 
focal plane has the important advantage of not requiring secondary mirror motion. 





For the purpose of this study we assumed that diversity defocus was introduced by 
translating the secondary mirror. This was simulated in the MTB setup by utilizing a 
fixed-lens WFS technique as shown in Figure 4. The lenses were designed and installed 
into a filter wheel with the defocus values indicated in Figure 4. 



Filter Wheel Locations & Defocus : 

1. +8.4 HgNg waves (note: Zernike basis) 


2 +6 4 waves 

3. -6.1 waves 

4. -10.6 waves 


Figure 4. Fixed- Lens WFS Setup for the MTB. 

WFS accuracy has been characterized earlier ~ using this technique and this is also the 
method of introducing diversity defocus using the JWST NIRCam. The technique also 
has the advantage simplifying the experimental procedure given that translating the InSb 
camera (to achieve defocus) is complicated by the bulky nature of its on-board dewar. 



Figure 5. Diversity Defocus Images with and without the MSA light shield. 


Sample mages collected at the focal plane of Figure 3 are shown in Figure 5. The top 
row shows the defocused images formed when the MSA is in place which clearly shows 
the MSA diffraction structure in the images. The bottom row shows the same defocus 
setting without the MSA. The sign convention for defocus is indicated at the bottom of 
Figure 5. 

As shown by the top row of Figure 5. the microshutters combine with the PSF complex 
amplitude and produce diffraction in the Fourier domain of the system stop. Since this 
image is relayed following the MSA. an additional pupil is formed and therefore a 3 ra 
direction plane exists before the image is formed on the CCD. Therefore, a multi-plane 





diffraction model correctly describes how the images are formed at the detector. This 
multi-plane diffraction propagation is illustrated graphically in Figure 6. 
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Figure 6. Multi-plane diffraction propagation. 


The calculation proceeds by first calculating the complex amplitude of the point spread 
function (PSF) in a defocused image plane at the microshutter: 

P s fca = 3{ pupil ca } = 3{ Ae' kn } , ( 1 ) 

where A is a mask defined by the telescope pupil geometry and (p d is a defocus phase 
term (the values used are listed in Figure 4). For JWST the mask A has an 18 segment 
hexagonal structure (see for instance Ref. 1- ) but for simplicity in our modeling we 
assume an un-obscured circular aperture. To properly incorporate diffraction by the 
microshutter light shield we next multiply the PSF complex amplitude by the 
microshutter light shield mask 

P s fca ~ P s fca x light shield mask . (2) 



Figure 7. Microshutter Array Light shield. 

The light shield mask is defined by an array of zeros and ones (1 = transmission; 0 = 
blocked) as indicated in Figure 7 which shows a microscope image of the light shield as 
well as its dimensions. 


The microshutters are assumed open over the region defined by the defocused image 
plane. If not, then the defocused images become severely under-sampled by the shutter- 
blocked portion of the images. An example of this is shown in Figure 8 which clearly 
indicates the consequence of not having all the micro shutters open. 







Figure 8. Diversity Defocus Image with un-opened Microshutters. 

After propagation to the MSA the beam is then re-collimated and stopped down by A, as 

shown in Figure 6. This calculation is modeled as an inverse Fourier transform to the 
complex pupil plane: 

pupil'ca =3“'{ psf' a } X A ; , (3) 

and then the complex amplitude at the detector is given by: 

P S f ca. final =%{P U P il 'j- ( 4 ) 

The measured intensity is the complex conjugate square of Equation (4). This calculation 
describes the multi-plane forward diffraction model illustrated in Figure 6 which is based 
on the MTB design of Figure 2. 

As an application of this approach consider Figure 9 which shows a direct comparison of 
the data with the forward diffraction model discussed above. The top row is the data 
(duplicated from Figure 5) and the bottom row are the corresponding calculated images. 
Clearly, the multi-plane forward model captures the physical effect of diffraction as seen 
by the diffraction blur the light shield in the of the bottom row of Figure 9. 



Figure 9. Diversity Defocus Image Data (top row) and Calculated Model Images 







The F /# variation with defocus is due to the non-telecentric nature of the experimental 

1 *3 

setup. The F/#’s were determined by a fitting procedure which matches the Fresnel 
zone structure of the images (concentric circles) with the specified amount of defocus 
shown in Figure 4. Initially, the F/#’s or defocus are not assumed to be known exactly 
but are used as starting values to the fitting procedure. The values shown in Figure 9 and 
Figure 4 are the final values determined by the fitting procedure. 


4. Image-based WFS using MOS Data 

One goal in this study has been to compare WFS results obtained from each of the data 
sets shown in Figure 5 to assess the impact on WFS accuracy due to the presence of the 
MS light shield. A multi-plane diffraction calculation was also examined to model the 
effect of the MS light shield on the intensity data in a defocused plane. Comparison of 
the modeled data to the testbed data shows good qualitative agreement as discussed in 
Section 3. In this Section we motivate the direct application of the image-based WFS 
algorithm to the data. Therefore, one considers a simple overlay of the MS light shield to 
the defocused intensity data, i.e., a direct multiplication of the light shield with the PSF 
intensity data. Therefore, no diffraction effect of the MS is considered since Equation (2) 
becomes: 


psf' 2 = psf 2 xlight shield mask , (5) 

as seen at the detector. A calculation comparison to the multi-plane model is made in 
Figure 10. 



data multi-plane simple overlay 


Figure 10. Diversity Defocus Image Data (top row) and Calculated Model Images 

The image on the far left is the data. The 2 nd image is based on the calculation described 
in Equations (1) to (4). The far right image is based on the calculation of only (1) and (5) 
. By comparison of the calculated images with the data above it is apparent that the 
diffraction effect captured by the multi-plane model is “blur” surrounding the MS light 
shield mask. However, at least to first order, the calculation described in (5) serves as a 
good approximation to the multi-plane approach. Physically, the reason this latter 
approach serves as a reasonable approximation is based on the observation that for a 
sufficiently defocused image, the image intensity begins to approximate the pupil. As a 
result, the effect of the MS light shield on the intensity data can be thought of as applying 
a binary aperture or obscuration at a pupil (the defocused image plane in this case 
approximates the pupil). As a result, one might expect that reasonably good WFS 
performance can be achieved without incorporating a multi-plane diffraction model into 
the forward Fourier propagations of the WFS algorithm. 





A number of image-based phase-retrieval techniques have been developed that can be 
classified into one of two general categories: (a) iterative-transform 14 13 or (b) 
parametric. !6, 1 Modifications to the original iterative transform approach have been 
based on the introduction of a defocus diversity function 16 ' ‘ 9 or on the input-output 
method 20 . Various implementations of the focus-diverse iterative-transform method have 
appeared in' 1 deviating slightly by utilizing a single wavelength and varying the 
placement and number of defocused image planes. Modifications to the parametric 
approach include minimizing alternative objective functions as well as implementing a 
variety of nonlinear optimization methods such as Levenburg-Marquardt, simplex, and 
quasi-Newton techniques." The current algorithm implementation strategy for JWST is 
discussed in these proceedings by Acton, et. al. 23 

In this paper, a further variant ~ 4 ' ~ 5 on the iterative transform method is applied that has 
the advantages of both the iterative-transform (high-spatial frequency) and parametric 
(high dynamic range) phase-recoveries. This is achieved by utilizing a "floating” or 
adaptive diversity function to incorporate feedback during phase-recovery to guide the 
PR process. A related approach has been discussed by Cohen, et. al.' 6 

The first case considered is the data set of the 2 nd row in Figure 5. i.e.. with no MS in 
place. The input data set to the WFS algorithm is thus 4 images - 2 on each side of best 
focus. Phase-retrieval results from this data are summarized by the GUI snapshot 
(graphical user interface) shown in Figure 1 1 . The results are then decomposed into a 

Zemike basis set with dominant components Z 6 =0.05 HeNe waves (astigmatism) and 
Z, 3 = -0.052 waves (spherical). The recovered phase is then used to recalculate one of 

the images of the data set (defocused image #4) to check consistency. This image is 
labeled as the MODEL image in Figure 11. X and Y cross-sections of the model image 
are also compared with X and Y cross-sections of the data. This latter comparison gives 
a qualitative consistency check on the assumed F /# and defocus values. These values can 
be fine-tuned using the GUI interface and the PR results subsequently updated. 
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Figure 1 1. Phase-Retrieval GUI and WFS Results without the Microshutter 
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Figure 12. Phase-Retrieval GUI and WFS Results with the Microshutter. 











The second case considered is the data set of the top row in Figure 5, i.e., with the MS in 
place. Phase-retrieval results from this data are summarized by the GUI image of Figure 

12. The results are again decomposed into the Zernike basis set. The MS data PR results 
differ slightly from the results of Figure 11 with components Z 6 =0.025 HeNe waves 

(astigmatism) and Z ]3 =-0.06 waves (spherical). The MS data PR results also have a 
larger component of first order 45 deg astigmatism. Z 4 . as well as second order 0-degree 
astigmatism, Z 14 . The recovered phase is again used to recalculate one of the images of 

the data set to check consistency. Although the PR results vary somewhat between the 
two cases, the results are still within general agreement given that the data sets were 
collected on two different days and the system could have been slightly disturbed while 
adding or removing the MSA. In summary for the low order PR difference between the 
two cases (with and without the MS) these two results are subtracted as shown in Figure 

13. The RMS difference between the two cases is = x/31 . 



Figure 13. Low-Order PR recovery difference with/without MS . 


An alternative comparison examines the PR recovery difference of the higher spatial 
frequency component by first subtracting the first 15 Zernike terms and then differencing 
the PR results with/without the MS. The phase difference between the two cases is 
shown in Figure 14. The RMS value is ~Xl 15. The structure of the higher spatial 
frequency component is better understood by calculating the PSD of the phase difference 
shown in Figure 14. This result is shown in Figure 15. The “cross-like” structure of the 
PSD indicates a remnant of the rectangular MS grid. Therefore, the higher spatial 
frequency component of the MS PR result contains the rectangular grid structure of the 
MS. As far as WFS and control is concerned, the rectangular component should be 
filtered to avoid any attempt by the control system to command this pattern onto the 
primary mirror. This filtering is effectively accomplished by decomposing the PR results 
into a lower order aberration basis set as indicated in Figure 12. 




Figure 14. High Spatial Frequency recovery difference with/without MS 


Figure 15. PSD of High Spatial Frequency recovery difference with/without MS 
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5. Summary and Discussion 

This paper examines PR results from a data set that includes a JWST MS array inserted 
into the defocused beam path. WFS results were reported on each of the data sets shown 
in Figure 5 to assess the impact on WFS accuracy due to the presence of the MS light 
shield. A multi-plane diffraction calculation was also examined to model the effect of the 
MS light shield on the intensity data in a defocused image plane. Comparison of the 
modeled data to the testbed data shows good qualitative agreement as discussed in 
Section 3. We have also shown that consistent PR results can be obtained by considering 
a direct application of the image-based WFS algorithm' 4 ' ' 5 to the data without any 
special modifications to the algorithm procedure. 

The MS element behaves as a diffracting element and introduces a rectangular artifact 
into the recovered phase when compared with PR results obtained from the 
corresponding data set without the MS. For the data sets considered in this study, the 
largest difference between the PR results with/without the MS array was due to the higher 
spatial frequency content. The RMS value of this rectangular artifact was approximately 
Z/ 15. By contrast, the lower-order PR results (represented by the first 15 Zemike 
coefficients) show that the PR difference is about a factor of two smaller or = /i/31 . 

A basic conclusion from the study is that image-based WFS using NIRSpec MOS data 
will produce PR results with two dominant components: (a) a low-order component, and 
(b) a higher spatial frequency component left as a remnant of the MS element. JWST 
will employ radius of curvature actuator control for each primary mirror segment and 
therefore the first component (a) will be the only component of the two which can affect 
the primary mirror figure through the WFS and control system. The higher spatial 
frequency component can be effectively filtered out by decomposing the PR results into a 
lower order aberration basis set as indicated in Figure 12. 

Additional work is underway which considers the effect of the JWST aperture on the 
diversity defocus data + the MS element. I.e.. one effect of the MS light shield is to 
slightly under-sample the data since some pixels are blocked in the data collection. The 
effect of this on segmented aperture recoveries is unknown. 

Another important effect not fully examined in this paper deals with the rather poor SNR 
of data collected by the MTB InSb camera. Many hundreds of frames must co-added to 
obtain narrow-band images of sufficient SNR to perform accurate WFS. For example. 
Figure 16 shows the InSb analog of the top row of Figure 5. The pixel size of the InSb 
camera is 30 jx square and therefore the MS light shield rectangular regions appear much 
smaller. Each image represents the sum of 100 broadband images after background 
subtraction and flat fielding. 




Figure 16. Broadband InSb Diversity Defocus Data Sets (A. G 1 to 5 //). 


Additionally, detector effects such as charge diffusion tend to be more pronounced using 
IR detectors and will likely further degrade WFS accuracy. Under such conditions the 
RMS differences between WFS results with/without the MS element could be 
substantially larger. The combination of blurring due to detector effects such as charge 
diffusion in addition to the inherent blur of broadband data will impact WFS accuracy28 
but generally only the higher spatial frequency component will be significantly impacted. 
Experiments are currently underway to obtain sufficient SNR in the IR image data 
utilizing a narrowband 2 n filter. Wavefront sensing error as a function of these factors 
will be reported in a future study. 
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